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Abstract: We construct uniform black-string solutions in Einstein-Gauss-Bonnet gravity 
for all dimensions d between five and ten and discuss their basic properties. Closed form 
solutions are found by taking the Gauss-Bonnet term as a perturbation from pure Einstein 
gravity. Nonperturbative solutions are constructed by solving numerically the equations of 
the model. The Gregory-Laflamme instability of the black strings is explored via linearized 
perturbation theory. Our results indicate that new qualitative features occur for d = 6, in 
which case stable configurations exist for large enough values of the Gauss-Bonnet coupling 
constant. For other dimensions, the black strings are dynamically unstable and have also a 
negative specific heat. We argue that this provides an explicit realization of the Gubser-Mitra 
conjecture, which links local dynamical and thermodynamic stability. Non-uniform black 
strings in Einstein-Gauss-Bonnet theory are also constructed in six spacetime dimensions. 
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1. Introduction 



In recent years it has been realized that, even at the classical level, the physics of black 
objects in higher dimensional spacetimes is much richer than in four dimensions. In order 
not to contradict observational evidence, the extra dimensions are usually supposed to be 
compactified on small scales. In such d-dimensional manifolds with n compact dimensions, 
black holes can either be localized in the compact dimensions, or the black hole horizon can 
wrap the compact dimensions completely [Q|. The static black holes which are localized in the 
compact dimensions have the horizon topology of a (d— 2)-sphere, S d ~ 2 , and are usually called 
caged black holes. For small values of the mass, these solutions look like the Schwarzschild 
black hole in d dimensions. In contrast, when the horizon wraps the compact dimensions, the 
horizon topology reflects the topology of the compact manifold. In the simplest n = 1 case, 
the single compact dimension is simply a circle, S . Then the second type of black holes have 
the horizon topology of a torus, S d ~ 3 x S 1 , and are referred to as black strings. 

In vacuum Einstein gravity, a uniform black string (UBS) is found by adding a trivial 
direction to a vacuum black hole in d— 1 dimensions. A static black string exists for all values 
of the mass, although it is unstable below a critical value, as shown by Gregory and Laflamme 
(GL) [0]. Following this discovery, a branch of non-uniform black string (NUBS) solutions 
(thus with a dependence on the extra dimension) was found perturbatively from the critical 
GL string in five [||, six [|j] and in higher, up to sixteen, dimensions ||. This non-uniform 
branch was numerically extended into the full nonlinear regime for d = 5 f| and d = 6 [Q, [| . 
In Ref. Q NUBSs were constructed for all dimensions up to eleven. 

All these results concern the case of Einstein gravity theory and its extension with abelian 
fields and scalars Q, [10]. However, for a spacetime dimension d > 4, the Einstein gravity 



presents a natural generalisation - the so called Lovelock theory, constructed from vielbein, 
the spin connection and their exterior derivatives without using the Hodge dual, such that 



the field equations are second order [11], [12]. Following the Ricci scalar, the next order term 
in the Lovelock hierarchy is the Gauss-Bonnet (GB) one, which contains quadratic powers 
of the curvature. As discussed in the literature, this term appears as the first curvature 
stringy correction to general relativity [13, 14], when assuming that the tension of a string is 
large as compared to the energy scale of other variables. The Einstein-Gauss-Bonnet (EGB) 
equations contain no higher derivatives of the metric tensor than second order and the model 
has proven to be free of ghost when expanding around flat space. Inclusion of a GB term in 



the gravity action leads to a variety of new features (see 15], [16] for recent reviews of the 
higher order gravity theories and their solutions). 

Although the generalization of the spherically symmetric Schwarzschild solution in EGB 
theory has been known for quite a long time [|l7|, |l8|], the issue of solutions with compact 
extra dimensions is basically unexplored. The only case discussed in a systematic way in the 
literature are the d = 5 UBSs in |ll| and their p— brane generalizations 1 in |p2|] . The results 



1 Note that these configurations are very different from the warped EGB solutions as discussed for instance 
in @, fjj. 
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there show the existence of a number of new features in this case, for example the occurrence 
of a minimal value of the black strings' mass for a given GB coupling constant a. The stability 
of these configurations was considered in fl23|| . However, the study of (asymptotically flat-) 
EGB black holes reveals that the case d = 5 is rather special [18] (see also Appendix A). 



Moreover, the results in j2q| for d > 5 UBSs with Anti-de Sitter asymptotics show that the 
properties of such solutions depend on the number of spacetime dimensions. Therefore it is 
of interest to study EGB black strings for d > 5 as well. 

The purpose of this work is to discuss the main features of the black string solutions in 
EGB theory for a number of dimensions between five and ten. Starting with configurations 
without a dependence on extra dimension, we first construct closed form solutions, by treat- 
ing the effects of the GB term as a perturbation around the Einstein gravity configurations. 
The nonperturbative UBSs are found by solving numerically the EGB equations with suit- 
able boundary conditions. For d = 5, the results in |T^] are recovered; however, no mass gap 
appears for d > 5 configurations. We shall address also the question of classical and thermo- 
dynamical stability of the UBSs. An interesting result here is the existence of black strings 
stabilized by the GB term in d = 6 spacetime dimensions. Black strings with a dependence 
on the extra dimension will be constructed for d = 6 only. 

The outline of the paper is as follows: in the next Section we explain the model and 
present general expressions for the physical quantities. The uniform black string solutions are 
discussed in Section 3 where we show approximate closed form solutions and results obtained 
by numerical calculations. In Section 4 we explore via linearized perturbation theory the GL 
instability of the black string solutions. Non-uniform black string solutions in d = 6 EGB 
theory are constructed in Section 5. We give our conclusions and remarks in the final Section. 
The Appendix contains a discussion of the basic properties of the EGB generalizations of the 
Schwarzschild solution in ci-dimensions, which helps us to better understand the properties 
of the black strings in d + 1 dimensions. 

2. The model 

2.1 The action and field equations 

We consider the (i-dimensional Einstein-Hilbert action supplemented by the GB term: 

' 1 ^(R+^L GB )d d x, (2.1) 



16ttG ./ V 4 

with 

L GB = R 2 - AR^R^ + R^ KT R^ KT , (2.2) 

where R is the Ricci scalar, g is the determinant of the metric, R a b is the Ricci tensor, R^ aKT 
is the Riemann tensor. The constant a in ( |2.l[) is the GB coefficient with dimension {length) 2 
and is positive in the string theory. 
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The variation of the action (2.1) with respect to the metric tensor results in the EGB 
equations 

1 a 

E^v = R^v — -^Rg^v + "jHfiv = , (2.3) 

where 

Hfj,v = ^{RpaKTRJ 7 ^ — 2R[ipvcrR PCT — 2RfiaR a u + RRpv) — -^LqbQiiv , (2.4) 

is the Lanczos (or the Gauss-Bonnet) tensor. 

A general feature of the EGB gravity is the presence of two branches of solutions, dis- 



tinguished by their behaviour as a — > [17]. However, in this paper we shall restrict to the 



branch of solutions which present a well defined Einstein gravity limit. 
2.2 Black strings in EGB theory: general formalism 

The simplest solution of EGB theory corresponds to the generalization of the Schwarzschild 
black hole. As we shall see, in (i-dimensions, such a configuration shares a number of basic 
features with the EGB black strings in d+ 1 dimensions. The basic properties of this solution 
are discussed in Appendix A. For example, it has an event horizon of S d ~ 2 topology and 
approaches asympotically the Minkowski spacetime background M . 

However, in this work we are interested 2 in solutions approaching asymptotically the d—1 
dimensional Minkowski-space times a circle, Ai^ 1 x S 1 . The line element of this background 
is 

ds 2 = -dt 2 + dr 2 + dz 2 + r 2 dQ 2 d _ 3 , (2.5) 

where the direction z is periodic with period L, r and t are the radial and time coordinates, 
respectively, while d^^_ 3 is the unit metric on S d ~ 3 . 

The physical quantities of a configuration that can be measured asymptotically far away 
in the transverse space are the mass M and the tension T in the direction of the circle. 
Similar to Einstein gravity, these quantities are defined in terms of two constants q, c z which 
enter the asymptotics of the metric functions 

5«^-l + ^4, ^-l + ^fc- (2-6) 
The mass and tension of a black string solution are given by 3 

M = ^§ [(d - 3)c t - c z ] , T =^[c t -(d- 3)c] , (2.7) 

where V^s is the area of the unit S d ~ 3 sphere. 

One can also define a relative tension n (also called the relative binding energy) 

TL ct-(d- 3)c z 
M (d- 3)c t - c z 



2 Note that, for small mass, the Schwarzschild- GB solution is expected to provide an approximation of the 
Kahiza-Klein caged black holes. 

3 For discussions of the computation of charges in EGB theory without a cosmological constant, see pq. 



- 4 - 



which measures how large the tension is relative to the mass. Uniform string solutions in 
vacuum Einstein gravity have c z = and thus a relative tension n = l/(d — 3). However, we 
shall see that this quantity is no longer constant when adding a GB term. 

The black string are higher dimensional black objects possessing an event horizon. For the 
solutions in this work, the event horizon is located at a constant value of the radial coordinate 
r = r h , with g u (rh) = 0. Similar to the black hole black string can be considered as a 

thermodynamical system, with entropy and temperature. The Hawking temperature of the 
solutions is given by 

T H = £, (2.9) 



with k the surface gravity. The general results in [27] shows that the entropy of a black object 



(i.e. also of a black string) in EGB theory can be written as a integral over the event horizon 

S=± j d d ~ 2 x^h(l + ^R), (2.10) 

where h is the determinant of the induced metric on the horizon and R is the event horizon 
curvature. 

The solutions should obey the first law of thermodynamics, which for Kaluza-Klein 
asymptotics contains an extra work term: 

dM = T H dS + TdL. (2.11) 

As we shall see, a Smarr relation relating asymptotic and event horizon quantities exists only 
for EGB black strings without a dependence on the z-coordinate. 

3. Uniform black string solution 

In this section we study EGB black strings which are translationally invariant along the 
z— direction. In Einstein gravity such configurations are constructed by adding a flat direction 
to any vacuum black hole in d — 1 dimensions. However, it is straightforward to check that 
this simple construction does not work in the presence of a GB term in the action, unless the 
solutions are conformally flat. Therefore the metric function g zz (and thus the proper length 
along the compact dimension) has a dependence on the radial coordinate r. 

3.1 The metric ansatz and field equations 

The uniform solutions are constructed in Schwarzschild-like coordinates, with the following 
metric ansatz: 

dr 2 

ds 2 = -b(r)dt 2 + J^+ 9(r)dfi 2 d _ 3 + a(r)dz 2 , (3.1) 
where a, b, /, g are functions of the radial coordinate r. 
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Substituting (3.1) into the field equations (|2.3D leads to a system of four equations for 
the functions a, 6, /, g. Fixing the arbitrariness of the metric gauge by taking g(r) = r 2 , these 
equations read 



a" + (?T + (d-3)-)^ 
2 r f 



-,'2 



«(/-!) 



+ (d-3)(d-4)- 
2a 7"v 



+ — 3)- 



a (d-3)(d-4)a 
4 r 2 / 



2(/-l)/-+2(d-5)(/-l)/- 
a ra 



. a 



+ (3/ - 1)/'- + 2(d - 5)(/ - 1)^ - (/ - 1)/^ + (d - 5)(d - 6) 



(/-I) 5 



(3.2) 



b" + ( h - + (d- 3)-)^ - * + (d - 3)(d - 4)^7^ + (d - 3)^ (3.3) 



a(--3)(d-4)6^ ?+ ^ 



4 r 2 f \ b rb 

+ (3/ - l)fj + 2(d - 5)(/ - 1)£ - (/ - 1)/^ + (d - 5,(<i - 6)^^ ) = 0, 

(d _ 2 ,4_^ + i ( ^4 ) + (rf _ 1)(ci _4)^ (3.4, 
rj zoo r a o r z j 

+ (d- 4)| ( - ^"^"^ ((d - 6)(d + 1)(/ - 1) + 2dr/') - + £)(6(d - 5)(/ - 1) 

+ r 2 /^ + 2r(3/ -!)£) + + + ^ (^/' + (<* " 5)(3/ - 1)) 



a& / r 2 a 2 6 2 r 2 a6 

2a" b 1 2b" a' \ 

+ 2 -^(rf b - - 2(/ - 1)) +%-Arf- - 2(f - 1) ) = 0, 

r A a b ' r A b a 'J 



>-»)(--^-(--S)^4)-^ (3-5, 

+ =((-- W -^ + ( ,-^(^^ ) + W - I ^)-a 

However, only three of the above equations are independent and one can show that ( |3.5D is 
in fact a differential consequence of fl3.2|) -(j3~4l). 

The equations under consideration are invariant under the following rescaling: 

a — > X 2 a, r — > Xr . (3-6) 

A dimensionless relevant parameter can then be defined according to j3 = a/X 2 , where A is 
some length scale. It is natural and convenient to choose A as the horizon radius of the 
black string and thus to define 

P=^. (3.7) 
r h 
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3.2 Boundary conditions and asymptotic solutions 

We require that our solutions approach asymptotically the background metric ( |2.5| ), which 
imposes a(r) = b{r) = f(r) = 1 as r — > oo. A black string has also a regular horizon at r = r^, 
where /(r/J = = 0, and a(r/ l ) > 0, /'(r^) > 0, b'(r^) > 0. Restricting to nonextremal 

configurations, local solutions in the vicinity of the event horizon r = can be constructed 
by expanding them in terms of r — rh as 

a = a + ai(r - r h ) + O (r - r h ) 2 , 

b = bi{r - r h ) + b 2 (r - r h f + O (r - r h f , (3.8) 
/ = fi(r ~ r h ) + f 2 (r - r h f + O (r - r h f . 

All constant in the near event horizon expansion are expressed in terms of the positive pa- 
rameters ao,b± (the expressions of a\, b 2 and f 2 are quite complicated and not enlightening; 
thus we shall not present them here). More important is the parameter /i which is fixed by 



= (3 2 (3d - U)d 3 d 5 dj + l2l3d 2 d A ds + 8d 4 - VP 
H ~ 4f3r h (pd 3 dld 5 + 2d 4 d 2 ) ' { ' ' 



with 



P = ^d\d\d\d\ + %fd 2 d\d\d^ + ?,2p 2 d G d b d\d 2 2 + Upd 7 d A d l 2 + G4d%, 



(for shortness, we have used here dt = d — k). Note that in the five-dimensional case, the 
obvious condition P > implies the existence of a minimum horizon size for a given value of 
the GB coupling constant 4 

3 = P < I ( 3 - 10 ) 
r h z 



i.e. for a > 0, one finds > v 2a. Since, as discussed in Ref. [19], the horizon radius 
is decreasing monotonically with the mass of the solutions, the relation ( |3.10D shows the 
existence of a minimal value of the mass for a given GB coupling constant a (note that a 
similar property is found for d = 5 Schwarzschild-GB black hole, see Appendix A). For d > 5 
there is, as far as we can see, no upper bound on f3. However, exploring the second derivative 
of the function f(r) around the horizon reveals the existence of another (negative) critical 
value: 

A = - ( rf -3) 2 (d-4) ' (3 ' U) 

with f 2 diverging as (3 — > (3 C . Although the existence of solutions with /3 < (3 C cannot be 
excluded, we could not find these disconnected branches. Then we conjecture that, similar 
to d = 5 case [p3f], for fixed, there exists a lower bound for the GB coupling constant, 



4 This is a nonperturbative effect, which cannot be seen considering the GB term as a small deformation of 
the Einstein gravity. 
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a > 



2r 2 

- ^-3)(d-4) • Seen oppositely, there exists an upper bound on the horizon radius (i.e. 
for the mass) for a fixed negative value 5 of a. 

For large values of r, the functions a, b, f admit the following expansion in terms of two 
constants q, c z 



a(r) 
f(r) 



1 + 



+ 



c t c z 



1 + 



d-4 1 2r 2 ( d ~ 4 ) 



+ . . . 



b(r) 



1 



+ 



c t c z 



d-4 1 2r 2 ( d ~ 4 ) 



+ ..., 



(3.12) 



ct 



r.d-4 



+ 



d-4 



CtC z 



2{d - 3) 2r 2 ( rf " 4 ) 



+ .... 



One can see that the GB coupling constant does not appear in the first terms of this expansion. 
However, we have found that asymptotic corrections specific to the GB interaction occur at 
higher orders in 1/r. 

Obtaining solutions extrapolating from (|3.8| ) to ( |3.12j ) allows to determine the values of 
the parameters an,&i,Q and c z as functions of and a, which are the input data in the 
numerics. 



3.3 A Smarr relation and scaling properties 

The thermodynamic properties of these solutions are determined by the asymptotic charges 
-the mass and tension, and the quantities on the horizon -the temperature and the entropy. 

The Hawking temperature Tjj of a black string as computed by using the definition 
in terms of the surface gravity or by demanding the regularity on the Euclidean section is 



4tt 



(3.13) 



The entropy 5 of a uniform black string is computed by using the general relation ( 2.10|) , 



V d _ 3 Lr d h ~ 3 ^h) 



4G 



l + -/3(d-3)(d-4) 



(3.14) 



Also, by integrating the identities 
R \ + j (Hi + \L gb ) = -rsz 



LA 

ab dr 



Sf (~1 + Hd 



3)((d-4)(/ 



Rl + j(H z z + lL GB ) 



LA 

ab dr 



l r d-3 a /. 



1) + /^)),(3.15) 
f (_1 + i a(d _ 3 )((d - 4)(/ - 1) + ff ,(3.16) 



between the event horizon and infinity, together with the field equations ( |2.3| ), one can prove 
an unexpectedly simple Smarr- type formula 6 . This formula relates quantities defined at 

5 Note that string theory requires a > and solutions with negative a may be unphysical. Nevertheless, in 
the study of UBS we have considered this possibility for the sake of completness. 



Note that the relation (3.15) can be used to derive the expression (3.14) for the entropy of UBS within 
the Euclidean approach to quantum gravity In this case, one has to supplement the action (2.1) with the 
corresponding boundary terms for Einstein and Gauss-Bonnet gravity. The tree level action / is regularized 



by subtracting the contribution of the background (2.5). Then the entropy is computed from F = M — TrS, 
with the free energy F = ThI, and M given by (2.7) As expected, the entropy computed in this way agrees 
with (Ell). 
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infinity to quantities denned at the event horizon: 

M(l - n) = T H S , (3.17) 

which also provides a useful check of the accuracy of the numerical solutions 7 . 

The rescaling ([O]) affects the thermodynamical quantities in the following way: 

T H ^^-, X d - 3 S 
X 

M -> A d_4 M, T -> A d " 4 T, (3.19) 

while a — > X 2 a, — > Xr^. 

We close this part by remarking that all UBSs in this work have an alternative interpre- 
tation as bubble solutions. The bubbles are found by using the analytic continuation z — > iu, 
t — > ix in the general line element ( |3.1[ ): 



dv 2 

ds 2 = -a(r)du 2 + b{r)d X 2 + + r 2 dQ 2 d _ 3 , (3.20) 

(where x has a periodicity (3 X = 1/Tjj). The properties of the bubbles can be discussed by 
using similar methods. For example, their mass and tension are fixed by the tension and mass 
of the corresponding black strings. 

3.4 Explicit solutions 

To our knowledge, for the equations (|3.2|)-(pO|) have no general closed form solutions. 



However, one can analyse the properties of the EGB uniform strings by using a combination 
of analytical and numerical methods, which is enough for most purposes. 

3.4.1 Einstein gravity UBSs 

The system of equations ( |3.2| )- (p7^) admit several exact solutions in particular limits. We start 
with a = solutions, in which case the family of {d — l)-dimensional black hole solutions of 
Schwarzschild-Tangherlini can be uplifted in d-dimensions: 

/( r ) = 6(r) = l-(^)^ 4 = / (r) , o(r) = l, (3.21) 

r = Th corresponding to the event horizon radius. The mass, tension, Hawking temperature 
and entropy of the Einstein gravity UBSs can be written in terms of as 

(E) _ LV d _ 3 (d-3)r d - 4 m _ V d - 3 r d - 4 {E ) _ (d - 4) (E) _ LV^ 3 



7 One can verify that this relation is trivially satisfied by the UBSs in Einstein gravity. However, the a — 
string solutions have another Smarr relation ]29] 

T H S = M d - 3 _- n , (3.18) 

which holds also for configurations with a dependence on the extra dimension. Due to the presence of a new 
length scale set by a, we could not generalize this relation to EGB theory. 
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These black strings exist for all values of > and have a constant relative tension n = 
\/{d — 3). Other properties of these solutions are similar to those of the (d — l)-dimensional 
Schwarzschild black holes. For example, their entropy decreases with the temperature. Thus 
the static Einstein gravity UBSs are thermally unstable. 

3.4.2 The small a-solutions 

An exact solution of the system (|3.2|) - (|3.4[) can be found in the limit of small a by treating 
the EGB configurations as a perturbation around the Einstein gravity. Here we have found 
convenient to take an ansatz with: 

a(r) = 1 + crai(r) + a 2 a 2 (r) + . . . , b(r) = f (r) (l + a6 x (r) + a 2 b 2 {r) + ...), (3.23) 
fir) = / (r) (1 + ah (r) + a 2 f 2 (r) + ...). 

In the first order in a one arrives at the system of linear ordinary differential equations 

„ 2(rf-3)-(rf-2)(^- 4 , d-3 f (rf-3)(rf-4) , 
JO a i + Yr 1 ~T ktl + ^ n 

1 , r ,\2d-4 

+ ^4 (f) ( d " 2 )( rf " 3 )( d " 4 )( d - 5 ) = 0, 

„ (2(d-3) + (d-6)(^) d " 4 ) . (2(d-3)-(d-2)(^) d " 4 ) . / x 

M' + — ~\ — % + — '—^ ^ -/ 3.24 

zr zr 

+ (d " 3 y~ 4) /i - (d - 2)(d - 3) 2 (d - 4)-L(^) 2 ^ 4 ) = 0, 

(d _ 4 + (d - 2)/ ) al + 2(d - 3)M + ^^/i 

r 

+ (d - 2)(d - 3)(d - 4)(d - 5)-^ (±f^) = o. 

zr r 

Similar equations are found for higher order functions, in which case they have a much more 
complicated form (and therefore are not shown here). When solving the equations for some 
order k in perturbation theory, there are four integration constants. These constants are 
choosen such that the corrected UBS metric still has an horizon at r = r^ and approaches 
the background ( pTq) asymptotically. 

Unfortunately, we could not find a general solution valid for any d and the equations has 
to be solved separately for any spacetime dimension. For example, in the first order in a, the 
solution reads 

2 1 2r h , , . 23 11 r h „ , , 1 7 5r h , 

oi(r) = H ^ + -4, Mr) = o tt, /lM = + o + — ^> ( 3 - 25 ) 

Sr^r 3r J 9r J 36r^r 36H 6r d 36rhr 36r^ 18r d 

for d = 5, and 

-M-St + lS 4l(r) = -l^-3. AW— 151 + ^,. (3-26) 
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in six spacetime dimensions. The expression of these functions becomes more complicated for 
d > 6. For example, the solution for d = 7 is 

4 /5\/3tt 15 6r 3 h lQy/3 ,2r + r h ^ 15 , . r h ,r h ^ 

«i( r ) = + ^ + 1- arctan(^^) - -j log(l + - + (-) 2 ), 



25 V r\ r 2 r$ r\ ^3r h 



r\ i r 



1 f _ VStt _ 3 ^(-138 + 20^ + 45 log 3) 63?i 
lrj "5/ (r)V r* r 2 + 10r3 + lOr^ 

+ ^ + 72-( 2 -— )( arctan (^^) + ^-log(l + 7 + ( 7 ) )) J, 
1 /(10^/3~7r- 78 + 451og3) 9rj| 33r£ 



5/ (r) V 10r 3 2r 5 10r 8 

SV3r h( 2r + r h a/3 r ft ,^2^ 

-^^( arctan (^^) + — log(l + 7 + (y) )) J, 

while for d = 8 one finds 

^)^ + ^|-^Ml + e) 2 ), (3.28) 

Mr) = — ( ~ — + r ' (151 ° g2 " 47) + ^1 + ^1 + (1 - 3(^) 4 ) A log (i + 

/o(r)\ 6r 2+ 9r 4 + 9r 6 + 6r 10 + 1 1 r ' ' 6r 2 g{ + 1 r ' ' 

I / rg(-32 + 151og2) 5rj 17r 8 5r 2 r, 
/l( )_ MOV 9^ + 3^ + 9^-3^ l0g(1 + ( 7 } ' 

Similar expressions have been found for d = 9, 10, without being possible to identify a general 
pattern 8 . 

This leads to following expression of the relevant quantities exhibiting first order correc- 
tions in (3 = a/r 2 : 

M = M^)(1 + Cl /?), T = T (E \l-c 2 (3), T H = T ( H E \l-c 3 f3), S = S^ E \l + c 4 /3), (3.29) 

where q are dimension-dependent positive numbers, e.g.: 

II 25 29 , , r 
ci = c 3 = — , c 2 = — , c 4 = — , for d = 5, 

36 36 18 

Cl= 16' C2= 16' C3 = 32"' C4 = ^2"' ford = 6 ' ( 3 - 3 °) 

_ 138 + 5^ -45 log 3 _ 5>/37r-9(18 + 51og3) _ 438 + 5y/&r - 45 log 3 

Cl ~ 50 ' ° 2 ~ 50 ' C3 ~ 150 ' 

2(288 + 5\/37r-451og3 r , „ 
c 4 = — , for d = 7, 

47-151og2 43 + 151og2 182-151og2 5(92 - 151og2) c , 
ci = - , c 2 = , c 3 = , c 4 = , for d = 8, 



8 The presence of logr and arctan(/iir + /12) terms (with /ii,/i2 dimension dependent parameters) is a 
generic feature of the d > 6 solutions. 
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M^ E \ T^ E \ Tft and being the corresponding quantities in the Einstein gravity as given 
by O- 

One can see that, in the first order perturbation theory, the mass and entropy increase 
with (3 while Hawking temperature and tension decrease. 

However, the situation is different when considering higher order corrections. For exam- 
ple, in the third order one finds the following expressions for the relevant quantities 9 : 

M = + —p + i?v + 6822513 n t = mi - ^ + «-f + 7796449 a (3.3D 

v 36 P 600 P 76204800 H h y 36 600 76204800^ ; ' v 1 

T„ = T^a --b + ™Lp - 2650383 3*) s = s^a + ™b + + 1063256 V ) 

H v 36 H 32400 P 228614400 H h y 18 H 32400 H 38102400 P h 

for d = 5, while the corresponding d = 6 expressions are 

M = M«(l + ^ - ™f - T = T«(l - 21/3 + ^ 2 - ^V), (3.32) 

v 16 3072 307200 y v 16 3072 P 307200 y v ; 

H v 32 P 3072^ 2457600 F ' v 32 H 1024 M 819200 ; 



It is easy to verify that the above expressions verify the first law ( |2.11 ) and the Smarr relation 
( pT7| ) up to order O ((B) 4 . 



One should also note that, while for d = 5 the mass and entropy always increase with 
j3, the situation is different for d = 6. Also, the above results make clear that the relative 
tension n cannot be constant for EGB black strings, since in the third order one finds 

1 _ 25 a , _67_ o2 , 7796449 q3 
n - 36^ ^ 600^ 76204800^ fn _ j _ r to oo\ 

?M I 11 I 127 .3 2 I 6922513 q3\ 101 fl_0 ' ^-"J 
^K 1 36 600 76204800^ / 

and 



31 a i 8995 2 _ 837467 o3 
i y ~Z 3072 ^ 307200 ^ 
Qfl _|_ _ 025-/92 _ 223547 03 
J| . iT 16^ 3072^ 307200^ 



- 16^ 1 3072 y 307200^ f j _ fi /o o 4 A 

n— 0/1 , 17 o 125 09 223547 IOT (3.34; 



The comparison of the relations ( 3.25|) -( j3.30 ) with those in Appendix A shows that, for small 



a, one can indeed treat an EGB black string as an uplifted EGB black hole. For example, 
in the first order perturbation theory, one finds the same qualitative dependence on a, with 
different coefficients, however. 

3.4.3 The large a-limit. A fixed point analysis 

In the limit where a — > oo (which is equivalent to taking only the GB part in the action), 
particular families of solutions can be constructed with the form a = r q , b = r p , f = fo where 
fo,p,q are constants. In fact, the underlying system of equations can be expressed in terms 
of A = ra'/a, B = rb'/b and f(r) and can be considered as a dynamical system in A, B, f 



The form of the metric functions are rather similar to (3.25), (3.26) 
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and in the variable r = log(r). Such systems admit usually one or several fixed-points. For 
example, in the pure Einstein-case, the equations can be set is the form 



r/' = 2(d-4)(l-/)-/(r) (A + B), rB'f = -{d-4)B, (3.35) 
(AB + 2(d - 3)(A + B) - 2(d - 3)(d - 4)) / = 2(d - 3)(d - 4), 



and admit only one fixed point corresponding to the flat space-time ( |2 . 5| ) . The corresponding 
equations for the GB theory are more complicated and we shall not present them here. 

For d = 5, these equations have p = q = 0, /o = 1 as the only fix-point. It corresponds 
to Kaluza-Klein space-time fl2.5|) ; obviously, this solution holds for all d. However, for d > 5, 
several solutions are available. More specifically, for d = 6, we find two more solutions, 
q = —1, p = 2, / = and a(r) = 0, p = 2, / = 0. The second is not a fix-point but 
nevertheless a (singular) solution. Note that b(r) ~ r 2 in both cases. 

In the case d = 8, we find more possibilities: 

g = 0,p = -2 5 / = l q = q, p = -4, f = --, q = 0, p = -1, / = 1, 

Q--— V--- f-~- q---V---f~- (3 36) 

g « 0.415, p = -5 + \/l3 rj -1.39, / ps 0.86, q « -3.82, p = -5 - « -8.6, / rj -0.28. 

A similar list exist for d = 7 and we expect the same situation for d > 8. One crucial difference 



between the d = 6 and d = 8 cases is that all the solutions ( 3.36 ) have p < implying that 
b(r) is divergent for r — > 0. The Ricci scalar is also singular at the origin. 

One may expect that if global solutions (i.e. for r G [0, oo]) of the full EGB equations 
exist in the limit a — > oo, these solutions would approach one of the fixed points above for 
r — > 0, and the Kaluza-Klein space-time ( |2.5|) for r — > oo. 

3.5 Numerical solutions 

The perturbative solutions above are useful since may provide some hints about the properties 
of solutions. However, this approach have some clear limits (e.g. it fails to predict the /3- 
bounds (jjUOD , (|3Tl|) ). 



Therefore we rely on numerical methods to construct the solutions within a nonpertur- 
bative approach. In this work, we have integrated the system of coupled non linear ordi- 
nary differential equations with appropriate boundary conditions by using a standard solver 
[pQ] . This solver involves a Newton-Raphson method for boundary- value ordinary differential 
equations, equipped with an adaptive mesh selection procedure. Typical mesh sizes include 
10 2 — 10 3 points. The UBS solutions have a typical relative accuracy of 10~ 6 . 

3.5.1 Pattern of solutions 

On general grounds, one expects a black string to share some of the properties of the corre- 
sponding EGB black holes in d— 1 dimensions. Therefore, it follows that the cases d = 5 and 
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d = 6 are special, since the GB term does not contribute in four dimensions, while the d = 5 
EGB black holes have distinct properties (see Appendix A). Moreover, one can see that some 
terms vanish in the equations (|3.2| )-(3,5) precisely for these values of d. 




The case d = 5 has been studied in [23]; our results in this case are in agreement with those 
of |23[] providing a good crosscheck of the numerical methods. The distinguished feature of the 
d = 5 is the existence of a finite range for the dimensionless parameter /3, with —1 < /3 < 1/2, 
with specific configurations approached at the limits of the /3-interval. However, the shape 
of the metric functions and the qualitative properties are similar to that found for d > 5. 
Therefore, for the rest of this part we shall focus on d > 5 solutions. 

Also, although most of the numerical data presented here corresponds to d = 6, 8, we 
have constructed solutions also for d = 7, 9, 10 (however, we did not study all these cases in a 
systematic way). Therefore we conjecture that UBSs in EGB theory exist for any d > 5 and 
we hope that our results catch the main features of the solutions available for an arbitrary d. 

For all configurations we have studied, a(r), b(r) and f(r) are smooth functions interpo- 
lating between the corresponding values at r = rh and the asymptotic values at infinity. A 
nonzero a leads to a deformation of the metric functions of a Einstein gravity black string at 
all scales and is particularly apparent on the function g zz = a(r) (note that this function may 
possess local extrema for some r > r^). As expected, for small values of /3, we have noticed a 
very good agreement of our results with the perturbative solutions ( 3.25 )-( ^28 ). The profiles 
of the metric functions of typical EGB black string solutions are presented on Figure |l| 

In practice, for a given value of d, we have constructed the full branch of UBS solutions 
both by setting a = 1 and varying the horizon radius and by taking = 1 and varying 
the GB coupling constant a. (Note that because of the scaling properties (p.6[), these two 
approaches are equivalent.) Also, for UBS solutions, the period L of the z-direction is an 
arbitrary positive constant and plays no role in our results. 
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Figure 2: The values of the metric function g zz an the horizon, a(r/ i ), and the derivatives of the metric 
functions at the horizon a'(rh),b'(rh), f'{i~h) are shown as functions of the dimensionless parameter 
j5 = a/ri for d = 6, 8 dimensions. 




As can be expected, when the value of the GB coupling is close to the critical (negative) 
value (3.11), the solver fails to provide reliable solutions. Also, for d > 5 we found no signal of 
a maximal value of a and of which may limit the domain of existence of black strings. For 
rh fixed and in the limit a — > (with both signs), the numerical solution smoothly approach 
the Einstein gravity UBS ( ^2l| ). 

From the point of view of numerics, the solutions can be characterized by the values of 
the event horizon parameters a(r/,), a'(r/i), /'(r^) and 6'(r^) (although only a(r/j) and b'(rh) 
are independent quantities). On Figure ||, these parameters are represented in function of 
the dimensionless parameter /3 for d = 6 (where = 0.5) and d = 8 (r^ = 0.2). It can be 
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Figure 4: The Hf component of the Lanczos tensor corresponding to a local 'effective energy density' 
is plotted for d — 6, 8 UBS solutions and several values of a. 



noticed that the dependence of these parameters on a is not always monotonic. Also, one 
can see that for a given (r^, a), a UBS solution with the proper asymptotics could be found 
for a single set of the horizon parameters a(r/ l ), b'(rh) only. 

Perhaps more relevant is the study of the branch of solutions obtained with varying 
and fixed a. The study of the limit — > needs a special attention since it cannot be 
approached by working with a finite (3. Our analysis of the black string in this limit reveals 
a completely different behaviour in the cases d = 6 and d = 8. As shown on Fig. || (left), for 
d = 6, the values a(r/j), a'(r/j), /'(r^) and 6'(r^) all tend to zero for — > 0. This suggest that 
the black string approach a configuration with a double zero of f(r), a(r) and b(r) at r = 0. 
The limit in the case d = 8 is easier to interpret: the metric approach the Kaluza-Klein 
background (2J>) on ]0,oo]. As shown on Fig. |3] (right) the quantities f'(rh),b'(rh) both 
tend to infinity for — > 0, indicating that the functions f(r),b(r) becomes discontinuous at 
the origin and that the flat metric (|2.5| ) is approached by the black string. Also, from the 
data in Figures 0, [|| one can notice the existence of solutions with g Z z( r h) > 1 and also with 
9zz{fh) < 1, for both d = 6 and d = 8. 

An explanation of the difference between the d = 6 and the d = 8 cases could be attempted 
by inspecting the corresponding fixed points and remembering that the study of the — > 
limit is equivalent to the study of the a — > oo limit, because of the scaling relations fl3.G|). 

In the case d = 6, a fixed point exists such that a = 0, / = 0, b = r 2 and the black string 
solution seems to be attracted by this fixed point in the region of the horizon for sufficiently 
small r/j. The limiting solution for = seem to extrapolate between the singular fix point 
for r — > and the background metric ([^) for r — > oo. This scenario is strongly suggested by 
Figure | (left). 

In the case d = 8, our classification of the fixed points ( 3.36| ) shows that b(r) becomes 
singular at the origin for all the solutions except for the flat spacetime background. As 
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Figure 5: The M, the relative tension n, the Hawking temperature Tjj and the entropy S of the 
d = 6, 8 uniform string solutions are shown as a function of the event horizon radius r/j, for a fixed 
value a = 1 of the GB coupling constant. For the data here and in Figure 8 we have set VasL/G = 1. 



a consequence, the 'Kaluza-Klein fixed point' appears as the most likely configuration into 
which the black string can be attracted for r e]0,oo] in the limit — >• 0. 

We close this part by noticing that, in the presence of curvature-squared terms, the 
modified Einstein equations leads to an effective stress tensor that involves the gravitational 
field G^ u = —aH^ u = T^ IU . Therefore, from some point of view, the quantity aH^(= —G^) 
corresponds to a local 'effective energy density'. However, this effective stress tensor, thought 
of as a kind of matter distribution, in principle may violate the weak energy condition. This 



property has been noticed in [31] for d = 5 black string solutions, being employed there to 
explain some properties of the black rings in EGB theory. As one can see in Figure f|, this 
feature of the five dimensional solutions is also valid for d > 5 black strings with a > 0. 

3.5.2 Thermodynamical aspects 

The various thermodynamical quantities were extracted from the numerical output. Our 
results satisfy the Smarr relation ( 3.17| ) with a precision of 10 -3 . We have also integrated 



the first law of thermodynamic J2.ll ) to crosscheck the mass obtained in this way for fixed 



values of L with the mass (2.7) computed from the asymptotic quantities c± and c z . These 
two computations of the mass are in excellent agreement, which provides a strong test of the 
accuracy of our numerical results. 

We start by presenting here the various thermodynamical quantities as functions of the 
event horizon radius r^, for a fixed value of the GB coupling a. Figure |5| (left) presents 
the entropy, Hawking temperature, mass and relative tension for d = 6 solutions. Although 
the data there has a = 1, this is the generic picture, since the results for other finite (and 



nonvanishing) values of a can be reconstructed by using the scaling (3.19). 



It can be observed that the temperature reaches a maximum at a particular value of 
while the entropy increases monotonically with r^; this implies the existence for d = 6 of two 
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Figure 6: The M, the relative tension n, the Hawking temperature Tu and the entropy S of the 
d = 6, 8 uniform string solutions are shown in units of the Einstein gravity uniform string solution 
(denoted by Mo, no, Tq, and So) as functions of the dimensionless parameter j3 = a/ri. 



phases of EGB black strings, one thermodynamically stable, the other thermodynamically 
unstable. (We will return to this property in the next subsection.) We note also that the 
mass and tension do not vanish in the limit — > 0. 

Figure || (right) presents the same quantities for d = 8 and a = 1. In this case, there is 
only one thermodynamically unstable phase since the temperature is always decreasing with 
r/j, while the entropy is always increasing. For d = 8, the mass and tension approach zero as 
r/j — > 0, confirming our interpretation that the solution approaches the flat spacetime ( |2.5| ) 
in this limit. 

In Figure 6 we illustrate the dependence of various thermodynamical quantites on the 
parameter ft = a /r? . There we have kept fixed the value of the event horizon radius and thus 
that data show how the solutions depend on the GB coupling constant a. The quantities 
there are given in units of the corresponding Einstein gravity uniform string solutions. An 
interesting feature here is the different behaviour of the solutions' mass. For d = 8, this 
quantity increases almost linearly with f3, while for d = 6 it approaches a constant for large 
(3. 



4. Stability of the EGB uniform black strings 

As shown by Gregory and Laflamme, for a given circle size, the Einstein gravity uniform 
black strings below a critical mass are linearly unstable. That is, at the critical point, there 
is a static mode breaking the translational invariance on the circle. This instability can 
be understood heuristically by using an entropy argument: for a black string horizon with 
S d ~ 3 x S 1 topology, there exist a length L above which it becomes entropically favourable for 
the mass to localize in a black hole with spherical topology. This connection between classical 
and thermodynamical stability has been made more precise with a conjecture formulated 
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by Gubser and Mitra (GM) |33|]. Their conjecture is that for systems with translational 
symmetry and infinite extent, a GL instability arises precisely when the system becomes 
thermodynamical unstable. 

Here we propose to find how a GB term will affect these results. The case d = 5 has been 
considered in |23]| . The results there show that GL instability persists for all allowed values 
of the GB coupling constant. 

4.1 Perturbations of the uniform metric 

In addressing the question of stability of the black string solutions presented above, we have 
found most convenient to use the approach of Ref. || and to consider in particular a NUBS 
metric ansatz 

ds 2 = _ b{r)e 2A(r,z) dt 2 + e 2B(r,z) + a(j . )d /\ + ^C^z) ^ 

where the non-uniformity is set in through the functions A,B,C. The uniform solutions 
discussed in Section 3 are recovered for A = B = C = 0. 

Then, following Q we perform an expansion of the functions A, B, C in terms of a small 
parameter e and consider a Fourier series in the z coordinate. In the leading order, we assume: 

X(r, z) = eXi(r) cos(fcz) + O (e) 2 , (4.2) 

X denoting generically A, B,C and k being the critical wavenumber corresponding to a static 
perturbation : k = 2n/L. 



We then substitute the form (4.1) in the general EGB equations (|2.3j) and expand A,B,C 



according to (4.2). In the next step, we consider the linear order terms there in the infinites- 



imal parameter e. The system obtained in this way is the relevant one for addressing the 



stability problem. Similarly to the pure Einstein case (for both flat ||] and AdS [32] black 
strings), the t — r component of the EGB equations allow to eliminate the field B\ in favor 
of the other fields and to reduce the problem to a system of two differential equations for A\ 
and C\. These equations are on the form: 

A'{ = ^A t + V>2^'i + V> 3 Ci + ^C[, C'l = tpiAt + ip 2 A' x + ipsd + ifiCi, (4.3) 

where the expressions of tpi, ifi are quite involved, so we prefer not writing them here 10 (they 
depend on the metric functions a, b, f and their first derivatives). 

Since we want to describe perturbations, the functions A\, C\ should obey the following 
conditions asymptotically: 

Ai(oo) = 0, Ci(oo) = 0. (4.4) 

The requirement that the perturbations obeying the linearized equations should be regular at 
the horizon leads to specific relations to be satisfied by A\{r^), C\{r^) and their derivatives. 



3 We have verified that the Einstein gravity equations of |1 are recovered in the limit a 
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The final expressions are really long due to the GB terms and we do not write them explicitely. 
We just mention that they are of the form 



A[(r h ) = A(r h) d,a;aa,bi;Ai(r h ),Ci(r h )), (4.5) 
C[(r h ) = C(r h ,d,a;a ,b 1 ;A 1 (r h ),Ci(r h )), 

where aQ,b\ appear in the UBS expansion (|3.8| ). Finally, the arbitrary normalisation in the 
linear system is fixed by 

Cx(r h ) = 1. (4.6) 

This leads to five boundary conditions which can be accomodated only if the parameter k 2 
assumes very specific values. The system of equations under consideration is therefore similar 
to a spectral problem. The sign of k 2 determines the stability: if k 2 is positive, there exists 
a GL unstable regime || for values of the wavenumber, say k, smaller than the critical wave 
number. If k 2 is negative, the system is stable. 



4.2 Numerical procedure and results 



The (linear) stability equations have been solved by using both the solver COLSYS 3Cj 
(which employs a Newton-Raphson method) and a two parameters shooting method. This 
has allowed to cross-check the results and to obtain the pattern of stability as accurately as 
possible. The shooting parameters can be chosen as the critical wavenumber k 2 and the value 
of the function A\ at the horizon. Together with the regularity conditions fl4.5| ), this leads to 
a Cauchy problem, suitable for a shooting method. 

The case d = 5 has been considered in p3| . The value of the dimensionless GB coupling 
constant j3 ranges 11 here in the interval [—1,1/2]. The results in [^3| show that the critical 
wavenumber k 2 is always positive and remains finite for f3 approaching the maximal value 
f3 = 1/2, and that k 2 diverges to infinity in the limit (3 — > —1. Our numerical results seem 
to confirm this behaviour, although we have found difficult to approach the limits of the 
(3— interval (see Figure 0). However, one can conclude that all d = 5 EGB black strings 
possess a GL instability. 

A main difference between the d > 5 and the five dimensional case is that the solutions 
exist for large positive values of (3. Our numerical results show that the critical wavenumber 



k 2 goes to a constant for large values of (3, as conjectered in [23]. However, we have found 
that this constant is negative for d = 6. Therefore, the EGB black strings are stable for 
large values of (3 in six dimensions. In other words, for a given event horizon radius, one can 
stabilize an EGB black string by taking a large enough value of the GB coupling constant. 

In the cases d > 6, the square critical wavenumber takes strictly positive values, i.e. the 
solutions remain unstable for any value of a. These results are illustrated in Figure for 
d = 6, 7, 8 and we expect a similar behaviour to hold for d > 8. That plot also shows that 



11 Note that the definition of the rescaled Gauss-Bonnet coupling f3 in this paper is different than the one /3 
used in 0, with $ = (1 - j3)/2. 
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Figure 7: The critical wavenumber is shown as a function of the dimensionless parameter {3 = ajr\ 
for d = 5, 6, 7, 8. The critical wavenumber diverges for {3 = — 1/3 in d = 6 (resp. /3 = —1/6 for d = 7 
and j3 = —1/10 for d = 8) and asymptotes to a constant for large values of (3. For d = 5, one finds 



— 1 < (3 < 1/2, the complete picture there being presented in Rcf. 23 



the value k 2 ^ = lim^oo k 2 increases with the number of dimensions of space-time. Also, as 
seen in Figure |?], the critical length L c of the extra-direction shows a complicated dependence 
on a. However, for large enough a and d > 6, L c increases with the GB coupling constant, 
approaching asymptotically a constant value. 

When the minimal value of (3 ( |3.11| ) is approached, we note that k 2 increases and likely 
diverges to infinity (note, however that the numerics becomes more difficult as (3 — > (3 C ). This 
behaviour, which also holds for d = 5 (see the Figure 3 in [23] and the footnote 11), appears 
also for the d > 5 data exhibited in Figure ^. 



One should notice that the perturbation (4.2) can be generalized to higher order, by taking 



into account the non-linear back reaction. For Einstein gravity black strings, the study of 
the perturbative equations in second order revealed the appearance of a critical dimension, 
above which the perturbative non-uniform black strings are less massive than the marginally 
stable uniform black string Q. It would therefore be interesting to solve the perturbative 
equations to second order also in the presence of a GB term, the only obstacle we can see at 
this moment being the tremendous complexity of the equations. 

4.3 Correlated stability conjecture 

We have seen that the family of EGB black strings present a dynamically stable phase in 
d = 6. However, for d > 6, we didn't find such a phase. In the following, we discuss this 
feature in relation with the thermodynamical stability of the EGB strings in d = 6 and d > 6. 




Figure 8: The entropy and the square critical wavenumber of d = 6, 8 EGB uniform black strins are 
shown as a function of the Hawking temperature. The d = 8 solutions have only one thermodynami- 
cally unstable phase. 



In the absence of matter fields or rotation, the study of the thermodynamical stability of a 
black string is related to the sign of the heat capacity 



C p = Th 



dS 
dfT 



H 



(4.7) 



Solutions with C p positive (respectively negative) are stable (respectively unstable). From 
the GM correlated stability conjecture [p5f| , it is expected that both a thermodynamically 
stable and unstable phases exist for d = 6 and that C p is always negative for d > 6. 

We have focused the analysis on d = 6 and d = 8, hoping that the dimension d = 8 
catches the main features for other values of d > 6. 

Figure || (left) presents the entropy S and the square critical wavenumber k 2 as functions 
of the Hawking temperature for d = 6. One can see that the solutions with k 2 < have also 
a positive specific heat (this branch has the index 1 in that plot), while the entropy decreases 
with Th along the branch with k 2 > (index 2 in the Figure ||). In agreement with the GM 
conjecture, the critical temperature where k 2 crosses zero coincides (within the numerical 
accuracy) with the temperature at the turning point of the UBS branches in a (S — Tjj) 
diagram. Therefore the correlated stability conjecture holds for the family of solutions under 
consideration. 

For higher values of the number of dimensions, there is no thermodynamically stable 
phase, as shown by Figure || for d = 8. As a result, the critical wavenumber k 2 takes strictly 
positive values. 

The different picture one finds for d = 6 can be understood heuristically by taking 
into account that the black strings should inherit some of the features of the BHs in d — 1 
dimensions. However, as discussed in the Appendix A, the d = 5 Schwarzschild black hole 
in EGB theory has a positive specific heat for large enough j3. For d > 5, the spherically 
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symmetric black holes are thermodynamically unstable, and thus also the corresponding d > 6 
uniform black strings. 



5. d = 6 non-uniform EGB black string 

As realized in [Q] , |J , the GL classical instability of the UBSs implies the existence of a new 
branch of solutions. The non-uniform black strings are known in Einstein gravity for all 
dimensions from five to eleven [|| , Q , ]?]] . These configurations are non- uniformly distributed 
along the z-direction and are continuously connected to the uniform strings, possesing the 
same horizon topology. However, for sufficient non-uniformity, the local geometry about the 
minimal horizon sphere (the waist) tends to a cone metric. The numerical results presented 
in [p5[| , Q suggest that, at least for d = 5,6, the non-uniform string branch merges with the 
caged black hole branch at a topology changing transition. 

Although all these solutions should present generalizations with a GB term, we restrict 
our study to the simplest case d = 6. This is also the most studied in Einstein gravity 12 , the 
numeric being easier for this spacetime dimension. Here we do not aim on a systematic study 
of these solutions, which would be a difficult taks. Instead, we want to establish the existence 
of EGB non-uniform solutions and to find their basic properties. Moreover, we have restricted 
our analysis to solutions with a positive value of the GB coupling constant (although we have 
verified the existence of solutions with a < as well). 

5.1 The ansatz and boundary conditions 

The non-uniform black string solutions presented in this paper are found within a metric 
ansatz employed also for the study of d = 6 NUBS in pure Einstein gravity 13 : 

r + r 

with ro a positive constant. The event horizon resides at r = where gu = and A,B,C 
stay finite. 

The equations for the functions A, B and C are found by using a suitable combination 
of the EGB equations, E\ = 0, E r r + E z z = and Ep = (with (p an angle on f^) which 
diagonalizes the Einstein tensor w.r.t. V 2 A, V 2 B, V 2 C (where V 2 = d rr + d zz ). Due to the 
GB contribution, these equations are much more complicated than in the case of Einstein 
gravity (with around 200 terms each equation) and we shall not present them here. 

The remaining equations El = 0, E r r — E z z = yield two constraints. Following Q, we 
note that setting E\ = E% = E r r + E z z = in V^E^ = and V^E» Z = 0, we obtain the 
Cauchy-Riemann relations 

d r Vi + d z V 2 = 0, d z Vi-d r P 2 =0, (5.2) 



12 The d = 6 NUBS have been constructed independently, by several authors J4|, |J, Jjj], with rather different 
methods. 



13 Although the metric ansatz (4.1) is useful in the perturbation theory, we could not employ it for the study 
of non-uniform solutions. 
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with V\ = ^—gE r z , V2 = ^—g(E^ — E")/2. Therefore the weighted constraints still satisfy 
Laplace equations, and the constraints are fulfilled, when one of them is satisfied on the 
boundary and the other at a single point ||. 

Utilizing the reflection symmetry of the non-uniform black strings w.r.t. z = L/2, the 
solutions are constructed subject to the following set of boundary conditions 

A\ =B\ =C\ = 0, 

lr=oo lr=oo lr=oo ' 

A\ n -B\ n = d ,d r A\ n = d r C\ =0, (5.3) 

\r=0 \r=0 u ' lr=0 lr=0 ' v ' 

® zA \z=0,L/2 = ® zB \z=0,L/2 = ^ zC \z=0,L/2 = °' 

where the constant cLq is related to the Hawking temperature of the solutions. Regularity 
further requires that the condition d r B\ r _ a = holds for the solutions. The boundary 
conditions guarantee, that the constraints are satisfied, since ^J—gE r z = everywhere on the 
boundary, and ^—g{E T r — E%) = on the horizon. 

One should note that, different from the Einstein gravity case, the functions A, B, C have 
a nontrivial shape for the uniform string, which has also do 7^ 0. 

The Hawking temperature and entropy of the black string solutions are given by 



TH = ^,S=^V 3 rU L dze^ 



1 + f£ (e- 2Co + e- 2B °r 2 (B ,*Co,* - 2C % - C , zz )) (5.4) 



where Aq(z), Bq(z), Cq{z) are the values of the metric functions on the event horizon r = tq. 



The mass and tension of a NUBS solution are still given by (2.7), where q and c z are the 
asymptotic coefficients in the metric functions gu = —e 2A ^ r ' z ^r 2 /(r 2 + 7q) and g zz = e 2B<yT,z \ 
respectively. 

A measure of the deformation of the solutions is given by the non-uniformity parameter 

a i 

^ \ 'Mnin / 

where TZ ma , x and 7Z m i n represent the maximum radius of the three-sphere on the horizon and 
the minimum radius, being the radius of the 'waist'. Thus the uniform black strings have 
A = 0, while the conjectured horizon topology changing transition should be approached for 
A -»• 00 |3§. 



5.2 Numerical methods 

While, as we shall see, their basic physical properties are similar to those of the Einstein 
gravity counterparts, the construction of NUBS in EGB theory is a more difficult task. This 
is due, perhaps, to the much more complicated form of the equations. 

To obtain these solutions, we solve the set of three coupled non-linear elliptic partial 
differential equations numerically, subject to the boundary conditions (|5.3j). The approach 
here is similar to that employed in |j] to construct d = 5, 6 NUBS in Einstein gravity. We 
employ a compactified radial coordinate r = r/(l + r) which maps spatial infinity to the finite 
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value r = 1. The numerical calculations are based on the Newton- Raphson method and are 
performed with help of the program FIDISOL/CADSOL |36|], which provides also an error 
estimate for each unknown function. For the solutions in this work, the typical numerical 
error for the functions is estimated to be lower than 10 -3 . 

The equations are discretized on a non-equidistant grid in f and z. Typical grids used 
have sizes 90 x 50, covering the integration region < f < 1 and < z < L/2. The horizon 
coordinate ro and the asymptotic length L of the compact direction enter the equations of 
motion as parameters. The results presented are mainly obtained with the parameter choice 
ro = 1, while the value of L is a-dependent. Also, due to numerical difficulties, our solutions 
have rather small values of the GB coupling constant a < 0.05. New features may appear for 
larger values, in particular for a/r 2 . close to one, a region which we could not explore with 
our methods. 

The most error-prone part of our numerical computation is to extract the asymptotic 
coefficients q, c z . Different from the Einstein gravity solutions, there is no Smarr relation 14 
which can be used to verify the global accuracy of the results. The existence of the Smarr 
relation ( |3.17 ) for UBSs was a consequence of the existence of a Killing vector d/dz, which 



implies the relation ( 3.1 6| ) . In fact, we have noticed that ( |3.17| ) does not hold for NUBSs, 



the deviation increasing with A. However, the Hawking temperature, the entropy and the 
non-uniformity parameter A are measured directly from the metric and thus from a numerical 
point of view they are 'cleaner' than the mass and tension. 

5.3 Properties of the solutions 

Similar to the case of Einstein gravity, a branch of non-uniform strings emerges smoothly 
from the uniform string branch at the GL critical point, where its stability changes. For a 
given value of a, keeping the horizon coordinate ro and the corrresponding asymptotic length 
L of the compact direction fixed, the solutions depend on the parameter do only 15 , which 
enters the Eq. ( |5,3[ ). Increasing this parameter, the non- uniform strings become more and 
more deformed, as quantified by the non-uniformity parameter A. 

In Figure |9] we exhibit the metric functions e 2A , e 2B , and e 2C for two typical solutions 
with different values of the GB coupling constant a and the same value of Tjj, ro and L. 
One can see that the shape of these functions is similar to that exhibited in the literature 
for a = 0. Moreover, one finds the same plots for the spatial embedding of the horizon into 
3-dimensional space. 

The A,B,C functions exhibit extrema on the symmetry axis z = at the horizon. As 
A increases, the extrema increase in height and become increasingly sharp. As A — > oo, one 
expects that the 'waist' of the NUBS would shrink to zero, locally forming a cone geometry, 
which should connect to the EGB caged black hole branch. Unfortunately, the numerical 
accuracy deteriorates for large values of A. 



14 This is not a surprise, given the presence of an extra-length scale fixed by a. We recall that no Smarr 
relation can be found even for the asymptotically flat black hole solutions in EGB theory. 



3 We recall that, within the ansatz (5.1), the UBSs have do 
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Figure 9: The metric functions e , e , and e of a typical d = 6 non- uniform string solution in 
Einstein-Gauss-Bonnet gravity are shown as functions of the coordinates r, z for a = 0.036 (left) and 
a = 0.01 (right). Both solutions have the same Hawking temperature Tjj = 0.125. 
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2.4 




Figure 10: The M, the relative tension n, the Hawking temperature Th and the entropy S of the 
d = 6 non-uniform string solutions are shown in units of the corresponding uniform string solution 
(denoted by Mo, no, To, and So) as functions of 1/(1 + A) (with A the non-uniformity parameter), for 
several values of the Gauss-Bonnet coupling constant a. 




Figure 11: Left: The M, the relative tension n and the entropy S of the d = 6 non-uniform string 
solutions are shown in units of the corresponding uniform string solution (denoted by Mq, no and So) 
as functions of the Gauss-Bonnet coupling constant a for a fixed value of the Hawking temperature. 
Right: The deformation parameter A is shown as a function of the Gauss-Bonnet coupling constant a 
for several fixed values of the Hawking temperature. 



Within our approach, we could find EGB configurations with reasonable accuracy 16 up 
to A ~ 3.7. 



We exhibit in Figure 1C the mass M, the relative tension n, the temperature Th and the 
entropy S of the d = 6 non-uniform string solutions, in units of the corresponding uniform 



Although solutions seem to exist for larger A, the numerical accuracy decreases drastically in that case. 
In particular, the constraint equations fail to be satisfied. 
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string solution, versus the non-uniformity parameter A, for several values of a. For any value 
we have considered, the key feature here is that the relative entropy grows with A, such that 
the NUBS are more entropic than the uniform critical string. Moreover, these solutions are 
also cooler than the critical GL strings. This is the picture one finds for a = 0. Also, one 
can see that as a increases, the derivation of the various quantities from the corresponding 
values of the uniform solution decreases. 

Turning now to the dependence of the solutions on a, we have constructed NUBSs in 
EGB theory starting with the d = 6 Einstein gravity solutions in || and slowly increasing 
the value of the GB coupling constant. Similar to the uniform string cases (Section 3) or 
to the black hole case (Appendix A), we have found numerical evidence that the known 
d = 6 Einstein gravity NUBSs also admit EGB generalizations. Keeping fixed the asymptotic 
length L of the compact direction, the horizon coordinate tq and the Hawking temperature 
Th, the solutions change smoothly with GB coupling constant. The results of the numerical 
integration in this case are shown in Figure 11 for several values of the Hawking temperature 
(thus of the boundary parameter <io). One can see that as the value of a increases, the 
non-uniformity parameter A descreases and the critical uniform solution with the given value 
of Th is approach. Along this branch, the entropy and mass decrease, while the relative 
tension increases. Therefore we observe again the tendency of the GB term to reduce the 
non-uniformity of solutions. It is tempting to link that to the occurrence of a negative local 
'effective energy density', ctH\ < 0, also for NUBS. This would imply the existence of an extra 
repealing force, which would reduce the string tension and will 'smooth' the configurations. 



6. Conclusions 

The phase structure of the d > 5 Kaluza-Klein theory with a single S 1 direction contains 
three 17 known phases of static solutions Q. These are the uniform black string phase, the 
non-uniform black string phase and the localized black hole phase. The main purpose of this 
work was to find how a Gauss-Bonnet term in the gravity action affects the properties of the 
first two phases of solutions mentioned above, with a S d ~ 3 x S 1 topology of the horizon. 

In the UBS case, we have answered this question by constructing the solutions using 
both analytical and numerical methods, for 5 < d < 10. Exact solutions were constructed 
by taking the GB term as a perturbation from pure Einstein gravity. The nonperturbative 
solutions were found by solving numerically the field equations. We argue that the presence 
of a GB term in the langrangian leads to some interesting new features, the cases d = 5, 6 
being special. For d = 5 there is a mass gap in the spectrum of black strings. In six 
dimensions, there are thermodynamically stable UBS solutions for some range of the GB 
coupling constant. Similar to the Einstein gravity case, the d = 5 and d > 6 black strings are 
thermodynamically unstable. 

17 However, at least for d = 5,6, in addition to these phases, there are also bubble-black hole sequences 
which we do not consider here. 
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We have also discussed the issue of Gregory-Laflamme instability for the UBS solutions 
with a GB term. All solutions were found to be unstable, except for the d = 6 case. There, 
for a given mass, the solutions with a large enough value of the GB coupling constant do not 
possess a Gregory-Laflamme instability. Thus our results provide also an explicit realization 
of the Gubser-Mitra conjecture, that correlates the dynamical and thermodynamical stability. 

However, one may argue that since the action (JjL^) is a one loop string theory approxima- 
tion, our results would be relevant for small value of a only. Therefore higher order corrections 
may change this picture, rendering the d = 6 solutions unstable as well. Moreover, for d > 8, 
higher order terms in the Lovelock hierarchy would enter the gravity action. Then we expect 
a more complicated picture in that case. 

One should also remark that the uniform solutions in this work can also be interpreted as 
black holes in a EGB-dilaton theory in d— 1 dimensions. The action of this model is found by 
doing a reduction with respect to the Killing vector d/dz, and has a very complicated form, 
with non-standard kinetic terms for the dilaton (see e.g. the Appendix A in Ref. pp[| ). (Note 
that the asymptotically flat five dimensional EGB-dilaton black holes will be thermally stable 
for some range of the parameters.) Following the usual procedure, one can add a U(l) charge 
to the d— 1 dimensional black hole solutions by boosting a black string in the (z, i)-plane and 
performing a Kaluza-Klein reduction with respect to the new z— direction. The properties 
of the resulting configurations are likely to be quite different as compared to the standard 
EGB-Maxwell black holes [41]. We hope to return on these aspects in future work. 



Concerning the NUBS branch of solutions, we have reported some partial results for the 
case d = 6 only. Our findings support the existence of EGB generalizations of the known 
Einstein gravity solutions. For the parameter range we have explored, we have notice a similar 
qualitative beaviour as in the absence of a GB term. In particular, our results suggest that 
the black hole and the black string branches merge at a topology changing transition also in 
this case, at least for small enough values of the GB coupling constant. There we have also 
noticed a tendency of the GB interaction to 'smoothen' the configurations, decreasing the 
deviations from the corresponding uniform solutions. 

To push forward our understanding of these issues, it would be interesting to construct 
the corresponding caged black holes. Given the results in this work, the dimensions d = 5, 6 
are of particular interest, and we expect the existence in these cases of caged black holes 
possessing a positive specific heat. 

Also, the EGB black string solutions should possess generalizations with matter fields, 
the case of Einstein-Maxwell-dilaton theory being perhaps the most important. However, the 
Harrison-type generation techniques used in to construct electrically charged solutions 
do not hold in the presence of a GB term, and one is constrained to approach this problem 
numerically. 
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A. The Schwarzschild black hole in EGB theory 

The generalizations of the Schwarzschild black hole solutions in EGB gravity were first discov- 



ered by Boulware and Deser [|17) , and Wheeler [|18|] , independently. The general solutions are 
classified into the plus and minus branches. In the a — > limit, the solutions in the 'minus 
branch' recover the ones in the general relativity, while there is no such limit in the 'plus 
branch'. Restricting to the 'minus branch' solutions, the line element of the Schwarzschild 
black hole in EGB gravity reads 

dr 2 

ds 2 = j^ )+ r 2 dtf d „ 2 -f(r)dt\ (A.l) 

where 

M = 1 + ^4)^4) (l - + f (d ~ 3)(d - 4) (l + (d - 3)(d - 4)^) i^- 1 ) . (A.2) 

In this relation, is the largest positive root of f(r), typically associated to the outer 
horizon 18 of a black hole, f(rh) = 0. The EGB Schwarzschild solution exists for all > 
and a > -2r 2 /(d - 3)(d - 4). As a -> 0, one recovers to leading order the Schwarzschild 
black hole expression, 

f(r) = 1 - Cff- 5 - (d "l ( f" 4) (^) d - 3 (l - a + ... (A.3) 

Also, as r — > one finds 



r ' Arf t r \ r 



Note that the usual form for /(r) in the literature is not given in terms of but rather in 
terms of a parameter [i = ?^ _5 (4r^ + (d — 3)(d — 4)a)/4. As seen bellow, this corresponds to 
fixing the mass of the solutions; however, the expression ( |A.2| ) facilitates a comparison with 
the numerical uniform black string solutions in which case the control parameter is the even 
horizon radius and not the mass, which is computed from the numerical output. 

In terms of the dimensionless coupling constant f3 = a/rf^, the Hawking temperature of 
a Schwarzschild- Tangerlini solution in EGB theory is given by 

= (rf-3) 4 + (d-4)(rf-5)/3 
H ATrr h 4 + 2(d-3)(d-4)/T 1 ' 



18 The global structure of this solutions is analyzed in |24]. Their generalizations with higher order curvature 



terms are discussed e.g. in 0], 
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with 



in the small f3 limit. Thus, for a > 0, the Hawking temperature of such systems appears to 
be suppressed relative to that of a Einstein gravity black hole of equal horizon area. 

At large distances (r 2 3> a) the EGB black hole behaves like the Schwarzschild solution. 
Once the event horizon radius is fixed (as is done in the numerical construction of the 
solutions), the parameter a enters the leading order term in the asymptotics 

/(r) = 1 _( 1 + 2 Wz^) (A „ 

which fixes the value of the mass: 

where V^_2 is the area of the unit S d ~ 2 sphere. One can see that the case d = 5 is special, 
given the existence of a mass gap M > 37ra/16G. The entropy of the Schwarzschild black 
hole in EGB theory, as computed from ( |2.10| ) is 

One can easily verify that the first law of thermodynamics dM = TndS also holds. 
This solution has a specific heat 

C =T ~ 2 >"~^ A + ( rf ~ 4 )( rf ~ 5)fl(2 + (<* - 3)(d - 4)/3) 2 

p H \dT H )~ 8(8 + 2(d-9)(d-4)f3 + (d-3)(d-4) 2 (d-5)f3 2 ) ' 1 ' J 

Without entering into details, we mention the existence of some substantial differences be- 
tween the thermodynamics of the d = 5 EGB Schwarzschild solutions and their Einstein 
gravity counterparts (here we shall restrict to the physical case a > 0). If the black holes 
are large enough, 3> y/a, then they behave like their Einstein gravity counterparts since 
C p < 0. A different picture is found for small values of r^, with Tjj — in that case and thus 
C p > 0. In fact, the specific heat changes its sign at = \fa. This implies the existence of 
a branch of small five dimensional EGB black holes which is thermodynamically stable (see 



the Ref. [15] for a review of these aspects). However, since C p < 0, the d > 5 solutions are 
thermodynamically unstable for all values of rh,ct > 0. 

The classical stability of the Schwarzschild black hole in EGB theory has been discussed 
in [^] for 5 < d < 11. The results there show that an instability occurs for d = 5, 6 at some 
large values of a (see also [0]). 
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